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Abstract 

We review recent developments in the theory and phenomenology of polarized 
structure functions. We summarize recent experimental data on the proton and 
deuteron structure function (71, and their impact on the understanding of polarized 
sum rules. Specifically, we discuss how accurate measurements of the singlet and 
nonsinglet first moment of gi test perturbative and nonperturbative QCD, and crit- 
ically examine the way these measurements are arrived at. We then discuss how 
the extraction of structure functions from the data can be improved by means of a 
theoretical analysis of their x and dependence, and how, conversely, experimen- 
tal information on this dependence can be used to pin down the polarized parton 
content of the nucleon. 
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1. The state of the art 



The state of the knowledge of polarized structure functions measured in deep-inelastic 
scattering has progressed during the last couple of years from the level of testing naive parton 
model ideas to the point of probing QCD at leading order (LO) and next-to-leading order 
(NLO). On the one hand, experiments at CERN and SLAC^~^) have accurately measured 
the structure function gi , which determines the leading twist contribution to the longitudinal 
polarization asymmetry, for proton and deuteron targets, over a reasonably wide range in x 
and at a pair of values of Q^. On the other hand, the perturbative behavior of gi in the (x, Q^) 
plane can now be quantitatively studied at NLO^-* , thanks to the recent determination of the 
full set of two- loop polarized anomalous dimensions^-'. This means that polarized structure 
function data can now be used to test perturbative QCD and obtain detailed information 
on the structure of polarized nucleons, while, conversely the data can be understood and 
analyzed using accurate NLO methods. 

Here we will review these recent developments, emphasizing the impact of the new data 
on the theoretical understanding of polarized structure functions. As it often happens, quan- 
tities that are easier to understand theoretically are the hardest to measure, and vice versa. 
We will take a theorist's viewpoint: we will start from the experimental results which have 
the most direct theoretical interpretation, and discuss them without, at first, questioning the 
procedure through which they have been obtained; we will then proceed to quantities whose 
interpretation is more subtle, and eventually discuss the procedure which is used to extract 
information from the data. 

In Sect. 2 we will review the determinations of the first moment of gi (as given by 
the experimental collaborations) and discuss in particular the behavior of the isotriplet first 
moment for which there is an absolute QCD prediction (the Bjorken sum rule). In Sect. 3 we 
will turn to the singlet component of the first moment, discuss its extraction from the data 
and its theoretical interpretation, which is complicated by the presence in this channel of the 
axial anomaly. In Sect. 4 we will then turn to the full x, Q"^ dependence of the structure 
function gi; we will discuss how it is described in perturbative QCD, and how it may be used 
to obtain information on the general structure of polarized parton distributions. Finally, we 
will come full circle, and show how this information can be used to improve the precision of 
the extraction from the data of quantities of simple and direct theoretical relevance, such as 
the moments of parton distributions. 

In this review we will concentrate on the physics of gi'^\ the discussion of other polarized 
structure functions being outside the scope of our treatment. It is however worth mention- 
ing that first measurements of (72, the other structure function which contributes to the 
deep-inelastic polarized cross section have been performed recently^^ While the theoretical 
interpretation of this structure function is worth studying for its own sake^-* , its contribution 
to the longitudinally polarized cross-section, on which we will concentrate, vanishes asymp- 
totically (as ^). In this contex, it is thus important mostly as a background: the recent 
experiments find that g2 is rather small, and essentially compatible with zero within errors, 
thus supporting the view that gi may be accurately determined by simply neglecting g2. 

2. The nonsinglet first moment and the Bjorken sum rule 

The main outcome of the recent precise experiments -^"^^ is a set of determinations of the 
first moment 
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Table 1: Summary of recent experimental determinations of Fi eq. (2.1) and ao eq. (2.4) (as presented 
by the experimental collaborations). The first error on Fi is statistical and the second systematic. 
All results hold at the average scale (Q^) (in GeV^) of the respective experiments. 



for proton and deuteron targets, displayed in table 1.* We will discuss later the theoretical 
input that goes into these determinations. If we take them at face value, however, they 
provide us directly with information that admits a simple theoretical interpretation, because 
Fi measures the nucleon matrix element of the axial current: 

M'als" = {t;p, s\iljaf.l5i^i\t;p, s), (2.2) 
ri(Q') = lj2^1C,^Q'^al = 1 [{e')Csmaim + Cmm<siQ')] , (2-3) 

where t indicates the target hadron (proton or deuteron, in our case) with mass M*, mo- 
mentum p'^ and spin s^, and the sum runs over all activated quark flavors; in the sequel 
when dropping the label t we will tacitly refer to a proton target. In the last step we have 
introduced the average quark charge (e^) = ^ 5^^=i , and the singlet and nonsinglet axial 
charges 



ao 



i=l 



(2.4) 
(2.5) 



i=l 



exploiting the fact that the Wilson coefficient functions C{Q'^) only distinguish between 
singlet and nonsinglet currents. 

The nonsinglet axial charge is particularly simple because the corresponding current is 
conserved (neglecting quark masses) and can therefore be decomposed into a sum of scale- 
independent contributions: 
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(2.6) 



^©(Q^ - Ql) (a„ + ad + as- 3ac) ± . . . 



* There exists also a direct determination of the neutron structure function and its first moment, obtained 
from scattering on a ^He target^°^ . We will not include this result in our discussion, because of its large statistical 
and systematic uncertainties — partly due to the problems related to scattering on a nuclear target '^^ — and also 
because a recent reanalysis of the same raw data has lead to a rather different determination of the corresponding 
structure function^ . 
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where (e^)^ is the average charge computed with Uf = k flavors, is the threshold for the 
k-th flavor, and all the scale dependence comes from the function ©(Q^ — <5|), which vanishes 
below threshold and is equal to 1 above threshold (having a nucleon target in mind in the 
last step we explicitly displayed thresholds for the charm and heavier quarks). 

The simplest nonsinglet scale-independent currents are then those related to light quarks, 
i.e. the isotriplet axial charge gA and the SU(3) octet charge ag 

9a^ CLu- Od; as ^ ttu + ad - 2as. (2.7) 

Using isospin symmetry, knowledge of Fi for proton and neutron targets determines the 
isotriplet combination: 

r{=i=rf-r? =CNs(Q')i « - aj) 

)-{^)-iS-l^f){^y (2.8) 

- (41.4399 - 7.6072n/ + iifnj) j + O(a^) -gA- 

The only assumption that goes into the derivation of eq. (2.8) is that of exact isospin symme- 
try, which implies = a^, = and = aJJ, where h indicates strange or heavy quark 
flavors. Furthermore, using isospin algebra again, gA in cq. (2.7) can be related to the matrix 
element of the axial current which induces beta decay, and thus to the beta decay constant 
of the same name, which is known rather accurately: gA = 1.2573 ± 0.0028. The relation 
eq. (2.8) between the beta decay constant and the isotriplet first moment of gi is known as 
the Bjorken sum rule^^^. In practice, the neutron structure function can be extracted from 
the deuteron one by using additivity of the deep-inelastic cross sections after correcting for 
nuclear effects^^^: 

rf^{l-1.5uo)^^, (2.9) 

where ~ 0.05 is the probability that the deuteron is in a D-wavc state. 

It thus looks as though the Bjorken sum rule (2.8) is an ideal place to test QCD in a 
clean and simple way: in principle, with a pair of determinations of Fj^^ at different scales 
both its normalization and its scale dependence may be measured, and thus isospin symmetry 
(which determines the value of gA) and perturbative QCD [which determines Cns(Q^)] may be 
simultaneously tested. Now, isospin symmetry violation in QCD is expected to be suppressed 
by powers of the current quark mass on a typical strong interaction scale; in particular in the 
unpolarized quark distributions it is at most of a few per cent ^^-^ ; the observation of larger 
isospin violation in this channel would presumably be the sign of nonperturbative physics 
and rather interesting per sc ^'^K 

On the other hand, Cns(Q^) is known up to order ct^, i.e., at the four-loop order, thereby 
allowing very precise tests of perturbation theory. In fact, at this order it might already be 
possible to see^^^ manifestations of the fact that the perturbative expansion in powers of 
as of QCD observables diverges^^^. A detailed treatment of this issue is beyond the scope 
of the present paper; its main implication for our purposes is that because of its divergent 
nature the perturbative series is ambiguous, in that different resummations of the series may 

differ by terms which are proportional to powers of . Because the values of physical 
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observables are unambiguous, these ambiguities must cancel against corresponding ones in 
other power-suppressed contributions, i.e. higher twist contributions. 
Indeed, in general eq. (2.8) should read 



(2.10) 



The coefficient Cht can be computed once a prescription for the treatment of the perturbative 
expansion is specified; the result for Fi is then unambiguous^^) . This program has never been 
carried through explicitly (there exists a proposal based on an explicit cutoff scheme which 
however would require a determination of Cht on the lattice). It is nevertheless possible to 
estimate the size of Cht phenomenologically, for instance using QCD sum rules; this leads 
to^o) Cht ~ -0.1 with an error of order 50% or more. Such a determination, however, is only 
meaningful if the ambiguity is substantially smaller than the value of Cht- Even though a 
qualitative argument"*^^) suggests that this is the case, an estimate-*^^) based on an analysis of 
the known terms in the expansion of Cns with the method of Fade approximants indicates 
instead that the ambiguity is of the same size as Cht- 

Be that as it may, these estimates give an indication of the level of accuracy at which 
perturbative QCD can be tested by comparing the Bjorken sum rule with the data. This 
comparison is shown in fig. 1, where F{^^ is displayed as a function of at various pertur- 
bative orders, with and without typical higher twist corrections, for different values of the 
strong coupling, and compared with an experimental value obtained"*^^) by averaging all the 
data of table 1 (plus some less precise earlier data) evolved to a common scale (we will discuss 
this evolution in the next section): 



r{=^(3 GeV^) = 0.164 ±0.011. (2.11) 



The comparison is striking: the result appears to be very sensitive to perturbative corrections 
— it would disagree with the theoretical prediction if this were computed at LO rather than 
NNLO — while being reasonably insensitive to higher twists and associated ambiguities. In 
fact, the value of as appears to be tested here to an accuracy comparable with the world 
average of all other determinations^^). In other words, if, rather than assuming isospin and 
using eq. (2.8) to measure ctg, we assign the value of ctg and test isospin, then the uncertainty 
due to the error on the value of OLg is already larger than the experimental uncertainty on the 
value of r(^^ quoted in eq. (2.11). 

This determination of r{^^, however, while being a theorist's dream, is a phenomenol- 
ogist's nightmare, and should be taken with extreme care. Firstly, the inclusion of nuclear 
effects in eq. (2.9) as a simple multiplicative correction factor may be an oversimplification: 
in a more detailed treatment of the deuteron wave function the correction is actually x de- 
pendent, and its effect on the first moment depends on the shape of the structure function 
g\{x)\ the effect has been argued to be non-negligible at the level of present-day accuracy, 
especially in the large-x region^^). Other nuclear effects, such as shadowing, which has a 
substantial effect"^^'^^) on the unpolarized counterpart of Ff^^ (the isotriplet first moment of 
F2), have been pointed out in this context^^\ but not studied systematically. 

Furthermore, it is clear that the reason why this test of perturbative QCD appears 
to be so sensitive is that the measurement is performed at a low scale, which exponentially 
magnifies the sensitivity to perturbative evolution. At such a low scale, however, other sorts of 
higher twist corrections, besides those discussed above, could be relevant, in particular those 
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Figure 1: Scale dependence of r{=^ eq. (2.10) [adapted from ref. ^^)]. The dashed curve corresponds 
to eq. (2.8) and the soHd curve to eq. (2.10) with Cut = —0.1; the dotted curve displays the 0{as) 
contribution to eq. (2.8), and the dotdashed is the asymptotic (a^ = 0) value. The three sets of 
curves correspond to the values as(M^) = 0.118 ± 0.007. The data point is as in eq. (2.11). 

related to the nucleon mass: even assuming that corrections due to transverse polarization 
(i.e. proportional to §2) are negligible, there are still kinematic higher twist corrections such 
as target mass corrections^^-*. Their effect on the first moment of r(^^ is^^-* comparable to 
the uncertainty in eq. (2.11) (which does not include it) at that scale. Moreover, such effects 
have never been studied systematically and could affect the extraction of all the first moments 
of table 1 from the data: even determinations with large (Q^) have several data points at 
low Q^. Finally, the result (2.11) is obtained by evolving the full first moment Fi, which in 
turn is obtained from data taken at various values of Q^. The corresponding uncertainty, 
related to perturbative evolution, could also be rather large at low scales: this will be the 
main subject of sect. 4. 

In conclusion, the fact that Vl^^ is only measured at relatively low Q'^ allows testing of 
the Bjorken sum rule only to about 10% accuracy, due to the large sensitivity to the value of 
as- This could be turned around: assuming the validity of the sum rule (i.e. of exact isospin) 
the value of can in principle be measured rather accurately. There is however a trade 
off, in that the determination of r{^^ at a low scale is affected by substantial systematic 
uncertainties, and the errors in table 1 may turn out to be over-optimistic, as we will see 
explicitly in Sect. 4. 

3. The singlet axial charge and the anomaly 

The singlet component of Fi can be extracted from the data in several distinct ways: by 
taking suitable linear combinations of different experiments, or by using SU(3) symmetry, or 
by exploiting the fact that the scale dependence of the singlet first moment differs from that 
of the nonsinglet. This scale dependence must also be taken into accoimt when comparing 
different experiments, and turns out to pose some interesting theoretical questions. We will 
discuss the phenomenological and theoretical aspects of these issues in turn. 
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3.1. The first moment and its scale dependence 



The simplest way of getting a value for the singlet first moment and the associated 
charge ao(<5^) [eq. (2.4)] is to subtract the nonsinglct contribution off each experimental 
determination of Fi. The axial charge is then found from eq. (2.3) by dividing out the singlet 
coefiicient function, which has been computed^^-* to NNLO: 



(3.1) 



The values listed in table 1 were obtained in this way. 

Neglecting all contributions from heavy quarks (we will come back to this assumption 
later) the nonsinglet is given by 

rr"°"=^CNs(g')(^A+^a8). (3.2) 

Just as qa can be obtained from the nucleon beta decay constants using isospin, as can be 
obtained from hyperon beta decay constants using SU(3) symmetry, with the result ag = 
0.579 ±0.025. This value is arrived at^^^ by a best fit based on the assumption of exact SU(3) 
symmetry; it could therefore by significantly affected by SU(3) symmetry breaking. This of 
course can only be introduced in a model-dependent way. One possibility is to introduce a 
parametrization of the hyperon beta decay constants which includes a term proportional to 
SU(3) breaking in the octet mass spectrum. A specific parametrization then leads to the 
value as = 0.40 ±0.19. An alternative option is to allow for current mixing effects, i.e. let the 
singlet current have a nonvanishing SU(3) nonsinglet matrix element, whose size can then be 
estimated within a model for the quark content of the octet baryons^-'^^ leading to ag ~ 0.4 
(and an uncertainty not smaller than the above) when SU(3) breaking is maximal. 

Even though the precise value of ag is thus affected by a large uncertainty, this has a 
rather small effect on the extraction of the singlet component, because its numerical coefficient 
is one third of that of the triplet and one fourth of the singlet. In fact, using these values 
of gA and ag in the decomposition (2.3) of the first moment into its singlet and nonsinglet 
components, it is easy to see that about 90% of the proton first moment comes from the 
isotriplet component, and of the remaining 10% about two-thirds are singlet, the rest being 
octet. For the deuteron the isotriplet contribution vanishes and the singlet and octet are 
partitioned as in the proton. This trivial numerology has some important consequences for 
the extraction of ao from the data: a) if the singlet is extracted from the proton data the 
uncertainty on ao will be about ten times larger than that on Fi ; b) an uncertainty of about 
50% in the knowledge of ag has the same effect as an uncertainty of a few per cent on the 
knowledge of the proton's Fi, i.e. it generates an uncertainty of about 10% in the proton or 
deuteron ao; c) the deuteron is in principle a better probe of the singlet component. This is 
reflected by the values of ao in table 1, derived neglecting SU(3) breaking. 

In order to meaningfully compare values of oq obtained in different experiments, they 
must be evolved to the same scale. Indeed, the singlet axial current is not conserved because 
of the axial anomaly, as we will discuss more extensively in the next section, so its matrix 
elements can acquire an anomalous dimension, which actually vanishes at LO, and has been 
determined up to NNLO^^^: 



ao(Q') = 



ao(oo). (3.3) 
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The full scale dependence of the singlet contribution to Fi [eq. (2.3)] is obtained by combining 
this with the singlet coefficient function eq. (3.1). 

It is then possible to determine qq from a global fit to the data of table 1 (rather than 
by simply averaging the determinations coming from each data set): if qa and ag are fixed 
using SU(3) symmetry as discussed above, ao(oo) is then the only free parameter. A global 
fit (including also as older data which however carry little weight due to their poor accuracy) 
with as{M^) = 0.118 ± 0.007 leads to32) 

ao(oo) = 0.29 ±0.04, (3.4) 

where the error is almost entirely statistical, because of the magnification of the statistical 
error on Fi when ao is extracted. 

Given the availability of independent determinations of Fi at different scales and for 
different targets, it is actually possible to relax the assumptions that went into the determi- 
nation eq. (3.4). First, it is possible to relax the assumption of SU(2) symmetry, and extract 
simultaneously ao and gA^'^^'- 

ao(oo) = 0.30 ±0.04; = 1-12 ± 0.10+°;^°. (3.5) 

The error on gA has been decomposed in statistical and systematic; the latter is entirely 
due to the uncertainty on ccg, the uncertainty on ag being essentially negligible. This shows, 
consistent with the discussion in the previous section and fig. 1, that the Bjorken sum rule 
can be tested at best to about 10% accuracy, due to the uncertainty on the value of ctg. 
Finally, it should in principle be possible to determine simultaneously gA, os and gq from 
the data: the proton-deuteron comparison fixes separately the isotriplet (i.e. gA) while the 
contributions of oq and as to the isosinglet could then be separated by their different scale 
dependence. This is however not yet feasible in practice because this scale dependence is only 
rather slight, while the sensitivity of the value of Fi to ag is weak; the latter fact, however, 
implies that even if we were to assume a 30% variation of ag due to SU(3) violation the value 
of ao (3.4)-(3.5) would hardly be affected. 

We now turn to polarized heavy quark contributions, which can be generated dynamically 
by assuming that the corresponding current matrix elements vanish on threshold. This then 
determines the various scale-independent nonsinglet contribution to eq. (2.6) in terms of the 
singlet at each threshold: for instance, the charm contribution is 

au + ad + as - Sac = [du ± ± a^] (Ql) = ao{Ql). (3.6) 

Thus, because of the small numerical value of oq, the error made neglecting heavy quark 
contributions to the nucleon Fi is indeed rather small. 

We can finally compare all the results of table 1 to one another and to the Bjorken 
sum rule prediction of eq. (2.8) by evolving to a common scale. The result, shown in fig. 2, 
demonstrates the mutual consistency of various experiments and the agreement with the 
prediction of exact isospin as tested by the Bjorken sum rule. 

3.2. The "spin crisis" and the anomaly 

The value of the singlet axial charge thus determined has attracted a good deal of 
attention because of its smallness. Small here refers to the expectation, based on the Zweig 
rule, that ao ~ ag, which the value in eq. (3.4) violates by several standard deviations. 
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Figure 2: Values of the first moment Fi eq. (2.f) from table 1 evolved to the common scale of 
= 5 GeV^. The solid lines correspond to the SMC experiments and the dashed lines to E143 
experiments. The dotted line is the prediction of the Bjorken sum rule eq. (2.8) (as the solid lines in 
fig. 1). 

However, as is scale- independent, while ao depends on Q'^ according to eq. (3.3): this raises 
immediately the question of the scale at which these two quantities should be compared. 
Indeed, the Zweig rule is a phenomenological rule of the naive quark model. Matrix elements 
of operators measured in hard processes should agree with the predictions of this model 
at some typical low hadronic scale: it is however clear that, at a sufficiently low scale, 
«o(Q^) [eq- (3.3)] grows as large as desired; in particular, it is of the same size as ag around 
0.5 GeV^. Whereas clearly perturbative evolution at such low scales cannot be trusted, 
nonperturbative estimates^^^ suggest that this scale dependence may actually be even stronger 
than perturbatively expected. 

More specifically, one may take the point of view that naive constituent quark observables 
are related to parton distributions measured in hard processes by a (generally x- and Q^- 
dependent) renormalization. A simple way of modelling this is to assume that the nucleon 
structure functions are given by combinations of constituent quark structure functions, each 
of which in turn is the convolution of a parton density inside the constituent quarks times 
the constituent quark density in the nucleon"^^'"^^-*. The constituent quark densities are then 
scale independent but target-dependent, while the parton densities satisfy the Altarelli-Parisi 
equations and should be universal. The model is tested by verifying this universality (for 
instance, it allows predicting the second moment of the pion parton distributions in terms 
of the proton ones) and seems to be in good agreement with experimental data^^-*. The 
violation of the Zweig rule, which is then by construction satisfied by constituent quarks, is 
attributed to their partonic substructure: the apparent contradiction with the quark model 
is thus removed, but the observed effect is not explained. 

A deeper understanding of the meaning of observed violation of the Zweig rule, and 
more in general of the first moment of gi, can be obtained in the QCD parton model. Matrix 
elements of operators are endowed with a partonic interpretation by identifying them with 
moments of parton distributions. The definition of the operators is however ambiguous since 
a renormalized operator can always be redefined by a finite renormalization: given a specific 
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definition of, say, the isotriplet axial current Js/^i, we can define a new current 

i5i=i = ^(««)i5j=i. (3-7) 

where Z[q;s(Q^)] = 1 + Z^^^asiQ"^) + . . .. The relation of the current to physical ob- 
servables will be modified accordingly: for example if r{^^ = C'ns(Q^) (is 7=1)5 then 
r{^^ = C^siQ'^) {Q'^) {j^'j^i) . This is, a priori, as good a definition of the axial cur- 
rent as the original one, and the ensuing ambiguity (factorization scheme ambiguity) only 
goes away at 00, or to all perturbative orders. It may, however, be fixed by physical 

requirements, such as the preservation of physical symmetries. In particular, if a current is 
conserved (as the nonsinglet axial current) it is protected against scale dependence: because 
dfij^ ~^ = then necessarily Q'^j^d^ji^ ~^ = 0, hence J5 (as a renormalized composite op- 
erator) does not depend on Q^. It follows that the redefinition of eq. (3.7) is not allowed if one 
imposes that chiral symmetry be preserved, because it would spoil this scale independence, 
and there exists a "natural" normalization of the current, namely, that which preserves chiral 
symmetry. 

The current matrix element can then be identified with the first moment of the polarized 
quark distribution: for instance the isotriplet charge is 

ai=i^gA^ [ dxAg/=i(x) = Aw(l,g2) _ Ad(l,g^), (3.8) 

where Au{x, Q^), Ad{x, Q^) are the polarized up and down distributions, whose A'"-th moment 
is given by 

Agi(iV,Q2)^ j dxx^-iAgi(x,g2); (3.9) 

and similarly for other nonsinglet charges. The conserved charge associated to is the 

total (isotriplet) quark helicity, so aj=i can be interpreted as the total (isotriplet) quark 
polarization, in agreement with naive partonic expectation. 

Things are more complicated in the singlet case; first, because there now are two parton 

distributions rather than one, i.e. the gluon distribution A (7 (a;, Q"^) and the singlet quark 
distribution 

AE(x,g2)^^Ag,(x,g2), (3.10) 

where Agj(a;, Q^) is the quark distribution of flavor i and Uf is the number of flavors activated 
at the scale Q^; and furthermore, because the current is no longer conserved, thus as a 
renormalized operator it depends on scale, and there is no unique "natural" way to normalize 
the quark singlet and gluon first moments, which will mix upon evolution. 

Indeed, the classical conservation of the singlet axial current is spoiled at the quantum 
level by the axial anomaly^^-*: 

^/^is = nf^ire^^'f"^F^,Fp,. (3.11) 

The singlet current is thus unprotected and requires an infinite renormalization, which induces 
a multiplicative scale dependence of the renormalized current of the form"^^) 

if, ^ l^{c^s)3l (3.12) 
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where t = ln((5^/A^). The anomalous dimension ^^{as) = li^^ct'i + • • • starts at two loops 
(because the nonconservation of the current is a one-loop effect) and depends on the specific 
choice of operator normalization; the scale dependence of gq induced by it was given in the 
MS scheme in eq. (3.3). 

Once this choice is made, however, the remaining normalization ambiguities are fixed 
by the anomaly equation. Indeed, eq. (3.11) holds as an equation between renormalized 
composite operators^^^ , and it therefore implies that the operator ire'^^^'^ F^j^Fp^ must mix 
with the current in such a way as to compensate (to all orders in perturbation theory) its 
scale dependence: 

n cv 

-Uf^tre'^^P'^F^^F,, = ^,{as)d^j^. (3.13) 

This condition fixes the normalization of the singlet gluon operator that mixes with the axial 
current, as well as the coefficient of this mixing — that there exists a scheme where eq. (3.13) 
is satisfied is the content of the Adler-Bardeen theorem"^^^ and modern versions thereof^^-*. 

However, there is still an obstacle in the identification of matrix elements of operators 
with partonic quantities: namely, we cannot simply identify the first moments AE(1, Q^) and 
Ag{l,Q'^) of the quark and gluon distribution [defined in analogy to eq. (3.9)] with matrix 
elements of two leading twist operators, because there is only one local leading twist (i.e. 
twist 2) operator with spin corresponding to the first moment of gi (i.e. spin 1). It follows 
that while it is natural to identify Ag{l, Q"^) (up to an overall constant, still arbitrary at this 
stage), with the quantity that mixes with the matrix element ao of according to eq. (3.13), 
i.e. 

-f^nf^Ag{l,Q^)^^,{a,)ao{Q'), (3.14) 

there is no unique way to identify ao{Q^) itself. Indeed, if the quark singlet and gluon 
anomalous dimensions are computed for higher moments, where a tower of quark and gluon 
operators (for all odd n > 3) are naturally identified with the corresponding moments of quark 
and gluon distributions"^^-* , and then analytically continued to the first moment, the result de- 
pends on the adopted regularization: if dimensional regularization is used throughout^^ , then 
the anomalous dimension of AE(1,(5^) coincides with 75(05); if instead infrared coUinear 
singularities are regulated by means of an explicit regulator^^) (such as putting the incoming 
particle off-shell) then the anomalous dimension of AS(1,(5^) = AS(1) vanishes. Note that 
in either case eq. (3.14) is satisfied, thereby fixing the constant of proportionality between 
Ag{l,Q'^) and the operator which mixes with oq . These two results correspond, respec- 
tively, to identifying the first moment of the quark distribution with the two eigenvectors 
of perturbative evolution: in the former case the singlet quark distribution is identified as 

AE(1, Q^) = ao(Q^), and in the latter case it is identified with the scale-invariant combination 
AE(l) = ao(Q2)+n/fjA^(l,Q2). 

This ambiguity can be pinned down only on the basis of physical requirements on parton 
distributions. Now, it may be shown'^'^^ that, in schemes where ao — S(l,(5^), soft con- 
tributions are partly included in the hard coefficient function, rather than being properly 
factorized in the parton distributions: this explains why these contributions are removed by 
an infrared regulator, which instead yields schemes where the quark is identified with the 
conserved eigenvector. Furthermore, in the latter schemes AS(1) can be directly shown^"*^) 
to coincide with the nuclcon matrix element of the canonical, conserved quark helicity oper- 
ator, the scale-dependence of ao then being due to a contribution to it from particle creation 
induced by the axial anomaly. In these schemes the Zweig rule expectation acquires a precise 
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meaning'^^\ since it predicts the approximate equality of two scale-independent quantities: 
ttg ~ AE(1). The (scale-dependent) deviation of ao from this value would then be explained 
by the gluon contribution to it^^^. Of course, whether this is actually the case rests with 
experiment: we will discuss this in sect. 4.2. 

The choice of factorization scheme in the definition of polarized quark and gluon distri- 
bution, besides affecting the physical interpretation of the singlet first moment of gi, has a 
significant effect in a finite-order computation of the scale dependence of Fi, and may thus 
substantially affect the extraction of this quantity from the data. Indeed, in the class of 
schemes where AS(1) is scale-independent, the singlet Fi is given by 

rf"- = l{e')CsiQ') (aS(1) - ^nfAg{l,Q')) (3.15) 

= ^ [C,(Q^)AS(1) + 2nfCgiQ^)Ag{l, Q^)] , (3.16) 
where the quark and gluon coefficient functions are given by 

C,{Q') = -^CiiQ') = -£^s(Q') (3.17) 

to all orders in perturbation theory. This is to be contrasted to schemes where AS(1, Q'^) = 
ao(Q2), soCg(Q2) = o. 

If the coefficient functions are determined at any finite perturbative order, eq. (3.15) 
implies that ao((5^) is not simply found, as in sect. 3.1, dividing F^'"^' by the coefficient 
function Cs((5^)- indeed, using the expression of Cq and Cg to order in eq. (3.16) one gets 

Cs(Q')<e^) rr^" - ^oiQ') = -Uf a':+'Ag{l, Q') + 0{a';+'), (3.18) 

where Cg^^ is the coefficient of ct^ in the perturbative expansion of CsiQ"^)- Thus all determi- 
nations of ao discussed so far are affected by an error of this size, and could only be improved 
by extracting AS(1) and Ag{l,Q'^) directly from the data (for example fitting eq. (3.16) to 
the data) and then computing ao{Q'^) — (AS(1) — ^2n/A(7(l, Q^)) . A more conservative 
point of view is that this is an intrinsic and unavoidable ambiguity of the computation. The 
size of the ambiguity of course depends crucially on the size of the gluon and decreases rapidly 
as the scale and the perturbative order increase. At order as, assuming that the quark re- 
spects the Zweig rule and the gluon makes up for the difference, and taking the value of ao 
from table 1, the correction at 3 GeV^ is around 15% of the value if ao ~ 0.5, and around 
50% if ao ~ 0.15. Note that this uncertainty is not included in the errors given in table 1. 

The scale dependence of ao induced by the axial anomaly is determined perturbatively 
by the anomalous dimension 75(0;^), which is a universal (i.e. target-independent) property 
of the axial current. It has been shown recently"^^-* that this universality actually persists 
beyond perturbation theory: the singlet axial charge satisfies a Goldberger-Treiman equation 
which relates it to the decay constant and irreducible coupling to the nucleon gf^^B of ^ 
singlet pseudoscalar state $ (defined in analogy to the tt coupling and decay constant in the 
usual Goldberger-Treiman relation): ao{Q^) = g^^B P^iQ^)- This pseudoscalar meson is 
however not a physical state of the theory due to the presence of the anomaly in this channel, 
so these couplings could only be measured on the lattice. However, the decay constant can be 
expressed in terms of the topological susceptibility, which is a property of the QCD vacuum. 



11 



Furthermore, all the scale dependence of ao comes from this quantity, and is thus contained 
in a universal, target-independent factor. It follows that if the explanation of the observed 
Zweig rule violation is in this scale dependence, i.e. in the smallness of the scale-dependent 
term the effect will be universal, and not specific of the nucleon. It is interesting to 
compare this possibility with an alternative proposal^^^ where the effect is instead related 
to a small value of the coupling g<^BB due to instantons. The effect would then be strongly 
target dependent, and have a different scale dependence. More data could thus shed light on 
our understanding of the QCD vacuum beyond perturbation theory. 

4. The structure function in the {x, Q^) plane 

We have seen in the previous sections that the scale dependence of Fi [eq. (2.3)] driven 
by the Wilson coefficient functions as well as by the anomalous dimension of the singlet axial 
current is quite large, that it substantially affects the extraction of current matrix elements 
from the data, and that it is actually responsible for a large part of the uncertainty on 
their determination. This suggests that a detailed understanding of the evolution of the 
full structure function gi in the (a;, Q^) plane is necessary in order to accurately assess this 
uncertainty, and pin it down as much as possible. Indeed, experimental data are taken 
within a limited range in x (0.003 < x < 0.7 for SMC and 0.03 < a; < 0.75 for E143) and at 
different scales in each x bin, with the lowest x points being taken at low and conversely 
(1.3 < Q2 < 48.7 GeV^ for SMC and 1.3 < < 9.2 GeV^ for E143). The determination of 
moments of gi requires thus both extrapolation in x and evolution to a common scale. The 
two problems are closely related because gi{x, Q^) is determined in terms of polarized parton 
distributions 

gi{x, = 1 [Cns ® AgNS + (e') (Cg ® AE + 2nfCg ® Ag)] , (4.1) 

(where ® denotes the usual convolution with respect to x). These evolve perturbatively 
according to the Altarelli-Parisi equations, so parton distributions at (xq, Qq) are causally 
determined from their values dX x > xq, Q < Qq. 

The values of the first moments in table 1 are arrived at by assuming the scattering 
asymmetry Ai to be scale-independent, i.e. by approximately evolving the data to a common 
scale on the assumption that the scale dependence of the polarized structure function gi is 
the same as that of the unpolarized one and extrapolating to small (and large) x by fitting 
a phenomenological shape to the last few data points. In order to critically examine these 
assumptions, we must study the small x behavior and scale dependence of gi in perturbative 
QCD. 

The small x behavior of polarized parton distributions 

Whereas the extrapolation of the structure function gi from the measured range to x = 1 
is under theoretical control and amounts to a small uncertainty, the extrapolation to small x 
is more subtle. This is due to the fact that gi must vanish identically for kinematic reasons 
at a; = 1; moreover, the form of its drop is already constrained by available data and agrees 
with expectations based on QCD counting rules On the contrary, the observed small x 
behavior has thwarted several times theoretical prejudice. In fact, at least in principle, the 
associated uncertainty is infinite, since the a; — > limit can never be attained (it corresponds 
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to infinite energy) and, because perturbative evolution proceeds from larger to smaller x, it 
is a priori impossible to exclude a growth in the unexplored small x region. 

The traditional expectation for the small x behavior of structure functions is embodied 
in Regge theory, which at some low scale should provide the small x behavior of parton distri- 
butions that are input to perturbative evolution . Dominance of Regge poles would predict 
that both the singlet and nonsinglet contributions to gi decrease or are at most constant as 
X ^ 0. Notice that in the unpolarized case Regge theory would have the singlet contribu- 
tion to F2 behave like a constant (or perhaps grow slightly, for a supercritical Pomeron) and 
the nonsinglet drop as \/a;; both the nonsinglet^^-' and singlet^^^ predictions are in excellent 
agreement with available experimental data. Because in LO F2 oc xq{x), while gi oc Aq{x), 
the Regge expectations for the singlet can be made consistent with the expectation^^) (partly 
based on QCD) that ^ ~ a;, if one assumes that both F2 and gi are constant or almost 

constant at small x at a low scale. The SMC data on g^, however, seem to indicate a rapid 
growth oi gi as X (see fig. 3a). One might have thought this to be due to a stronger 
increase of the singlet at small x, perhaps due to Regge cuts, but this is belied by the SMC 
data on gf, which should then look like the proton in the small x region, whereas in actual 
fact it is negative there (fig. 3b). 

Now, even if at some low scale parton distributions behave according to Regge expec- 
tations, they will be modified by perturbative evolution. In particular, it has been known 
for long that LO perturbative evolution eventually leads to a rise of gi at small x'^^'^'^^ of 
the same form as the rise of the unpolarized structure function ^2^^-*; the main difference is 
that whereas F2 is positive-definite gi is not, so the rise may correspond to gi growing either 
large and positive or large and negative. A closer look reveals both analogies and differences 
between the polarized and unpolarized cases^°) . The small x behavior of parton distributions 
is dominated by the rightmost singularity of perturbative anomalous dimensions in the space 
of moments. In the polarized case, this singularity is located at = for the nonsinglet as 
well as for all entries in the matrix of singlet anomalous dimensions, and it is already present 
at LO. In the unpolarized case it is located at A = 1 for the nonsinglet and at A = for the 
singlet; in the latter case the singularity is present at LO in the gluon anomalous dimensions, 
but only starts at NLO in the quark ones [defining the moment variable A according to 
eq. (3.9)]. 

As a consequence, both the singlet and nonsinglet polarized distributions grow according 

to^o) 

A/(x,Q2)^^e2^/^ (4.2) 

where a = ^/^, p = \/C/Cj C = ^5 C = ^"^q^] , and xq is a reference value of x such 
that the approximate small x form of the anomalous dimensions is applicable for x <xo and 
> Qo- nonsinglet case. A/ = AqfNS and 

In the singlet case the quantities that display this growth are the linear combinations of quark 
and gluon v"^ = AS -|- C^Ag, with 



^" = 2(1±^/1-|^ ; 7l = 7Ss 5±4^1-|^|. (4.4) 
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Notice that this is a scheme independent result. The behavior of the unpolarized nonsinglet 
is the same as that of the polarized singlet (but with different 7), while the behavior of the 
unpolarized singlet differs, first because the quantity which grows according to eq. (4.2) is 
xg{x, Q^) rather than Ag{x, Q^), and furthermore because the growth is driven by the gluon, 
i.e. the unpolarized quark is proportional to the unpolarized gluon (and of order as with 
respect to it). 

This implies that eventually the singlet distributions indeed dominate gi at small x. 
In both singlet eigenvectors the quark and gluon have opposite sign, and asymptotically 
AE > if AE(a:o, Qo) > —C~^ Ag{xo, Qo): hence, for most plausible starting quark and gluon 
AE grows large and negative and A^r large and positive. As a consequence, asymptotically 
gi at small x must be large and negative. The nonsinglet distribution however also grows 
asymptotically; the growth has the same form as for the singlet distributions, only with a 
smaller value of 7/, its sign is the same as that of the starting distribution. Therefore, the 
value of X where the singlet actually does dominate could in practice be extremely small. 
This is to be contrasted to the unpolarized case, where the growth of eq. (4.2) is down by a 
power of X in the nonsinglet compared to the singlet, so that the singlet already dominates 
around x ~ 10~^. 

Of course, the perturbatively generated growth of eq. (4.2) is only visible provided it 
dominates over a possible growth of the starting parton distributions: even though Regge 
theory does not favor such a growth, it cannot be firmly ruled out. For instance, if parton 
distributions rise as Af{x,Qo) ~ then the evolved distributions will behave according 
to eq. (4.2) only for p < but for smaller values of x they will behave as^°^ 



Af ^x 



(4.5) 



i.e. they reproduce the boundary condition up to an x-independent correction; in the singlet 
case if the starting quark and gluon have different small x-behavior the most singular one 
dominates asymptotically. In such case, the onset of the asymptotic small x behavior of 
eq. (4.5) will be yet slower^^^: firstly, because the small x solution to the evolution equations 
is then dominated by the singularity of the boundary condition rather than the singularity 
at AT = in the anomalous dimension, therefore, approximating the anomalous dimension 
with this singularity is less accurate; also, the two singlet eigenvectors only mix due to the 
x-independent correction in eq. (4.5), thus the dominance of the most singular behavior only 
sets in rather slowly. 

It is interesting to see how these LO predictions are modified at higher orders. At NLO 
all polarized anomalous dimensions have a singularity. This leads to a correction to the 
LO of the form^) 



A/^LO(^^q2) ^^^^ ^^(^)3 ^^^^^2) _ ^^(Q2))j Af^^{x,Q^), 



(4.6) 



where A/^*-'(a;, Q^) is given by eq. (4.2), and the coefficients e are explicitly 
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Figure 3: Plot of gi{x) for (a) proton and (b) deuteron. The crosses are SMC data and the diamonds E143 
data. The curves correspond to a fit^^Ho the proton data only with valence-like nonsinglet and (i) "maximal" 
gluon and steep input quark singlet or (ii) "minimal" gluon and flat input quark singlet. 



corresponding to a further rise proportional to ^^/^ of the parton distributions at small x, 
with coefficient cns/Tns ~ ~ |- Notice that the form of the small x eigenvectors is 

unaltered, and that these results are scheme-independent. 

At NLO the Mellin transforms of the quark singlet and nonsinglet and gluon coefficient 
functions also acquire a singularity at N = 0, which is however only of order and therefore 
does not contribute to the asymptotic small x behavior eq. (4.2). This singularity corresponds 
to a In i rise (with positive sign) of the functions C{x, Q^) in eq. (4.1). At small x, this leads 
to a rise of gi even when the various parton distributions are fiat. Even though this rise is 
weaker than that induced by perturbative evolution, it could dominate at low scale where 
evolution effects are small; its precise form is scheme-dependent, however. 

At yet higher orders, the expected generic behavior of anomalous dimensions is 7 ^ 

(ct^/AT^"'"^). The small N behavior of the NLO coefficient functions suggests that these 
might instead behave as C {as/N"^)^ ^, which would guarantee scheme independence 

of the coefficients of the singularities in the anomalous dimensions at least in the nonsinglet 
case. Indeed a resummation of these singularities in the nonsinglet anomalous dimensions 
has been proposed long ago^^^. In the absence of appropriate factorization theorems, it is 
unclear whether this resummation, which reproduces the known NLO leading singularity, 
correctly reproduces the leading singularities of the nonsinglet anomalous dimensions to all 
orders. If it does, it implies that the polarized and unpolarized nonsinglet quark distributions 



15 



will display the same small x-behavior, in contradiction to Regge theory. More generally, one 
would expect^^^ the summation of logarithmic effects in In ^ induced by these singularities 
to lead to a rise of parton distributions, which should be asymptotically power-like; this rise 
should set in significantly faster than in the unpolarized case, due to the stronger nature of the 
singularity, and thus, contrary to the unpolarized case^^\ it could set in before perturbation 
theory breaks down. 

The combination of all these effects makes it rather hard to infer the asymptotic small 
X behavior of polarized parton distributions from present-day data: a HERA-like kinematic 
coverage would be necessary to disentangle the various effects. The only firm prediction at 
this stage is that parton distributions at small x will certainly grow (positive or negative) at 
least as fast as eq. (4.2), including the nonsinglet distributions, and that therefore a simple 
phenomenological extrapolation based on Regge expectations is not adequate at least in 
principle. More detailed statements can be made by studying the structure of the data in the 
(x, (5^) plane. 

^.2. The evolution of gi{x,Q^) 

The only way of arriving at a precise and reliable determination of the moments of gi 
is to describe its evolution in the (x, Q^) plane in terms of the evolution of polarized parton 
distribution. The fact that these evolution effects are large not only implies that they must 
be included in the data analysis, but also that they might provide useful information on 
the polarized parton content of the nucleon, and specifically on its gluon content. In fact, 
one would expect gi to be a very sensitive probe of Ag — much more than the unpolarized 
structure function F2 is a probe of g — since, even though the coupling of the gluon to gi is 
formally NLO, thanks to the anomaly it does not decouple asymptotically^^-*. In comparison 
to the unpolarized case, however, the analysis is complicated by the fact that the singlet 
and nonsinglet display a similar perturbative behavior in the small x region: this implies 
that only combining proton and neutron (or deuteron) data is it possible to disentangle the 
contributions of the various parton distributions to gi, at least using data with the current 
relatively restricted kinematic coverage. 

Thanks to the recent determination^^ of the full matrix of two-loop anomalous dimensions, 
it is possible to study the evolution of parton distributions at NLO. The computations of 
ref. ^) are performed in the MS scheme, where the gluon decouples from the first moment of 
gi and thus is not properly factorized, as we discussed in sect. 3.2. However, since the quark 
and gluon coefficient functions are known"^^'^^-* in schemes where the gluon does contribute 
to gi, the corresponding anomalous dimensions may be determined^-* from those of ref. ^\ 
A determination of parton distributions can then be made by fitting to the data, provided 
one makes some simplifying assumptions on the form of the starting parton distributions.* 
Specifically, a simple polynomial parametrization of parton distributions may be adopted, 
which entails the assumption that the asymptotic small x behavior of parton distributions 
starts setting in at the values of x explored at present: 

A/(x, Ql) ^N{aj,/3j, af) ry/x"/ (1 - x)^f (1 + afx) (4.8) 

* For a review of LO polarized parametrizations see ref. . A NLO fit has been presented in ref. , but 

in a scheme where the gluon does not couple to gi. Approximate NLO fits (prior to the full determination of 
NLO anomalous dimensions) in properly factorized schemes are in refs. 58,50)_ ^ £^jj NLO computation in a 
variety of physical schemes is in ref. . 
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Figure 4: Plot of (i) gi{x) and (ii) xg-i{x) for (a) proton and (b) deuteron. The crosses are SMC data and the 
diamonds E143 data. The curves correspond to a NLO fit^Ho all data. 



{where N {a, P, a) is a normalization factor, fixed e.g. as N{a, P,a) J^dx x°'(l — x)f^(l + ax) = 
!)• 

Two fits^*^-* to the proton data only of parton distributions of this form are displayed in 
fig. 3, together with the corresponding predicted deuteron structure function. The nonsinglet 
is assumed to be valence-like (q;ns = +0.2); the gluon is assumed to be fiat {ag — 0) and the 
singlet quark is either steep (ctg = —0.5, fig. 3i) or fiat (ctg = 0, fig. 3ii). The starting gluon 
is fixed either by requiring AS(1) to satisfy the Zweig rule, so that the fit then forces a large 
gluon component to compensate ("maximal gluon", fig. 3i), or by symply taking the gluon 
distribution to vanish at the starting scale ("minimal gluon", fig 3ii). The quality of these 
two fits does not differ in a statistically significant way: while a larger gluon corresponds to 
stronger evolution effects, the proton data alone do not allow us to fix its size, or the small 
X behavior of the various parton distributions. Notice, however, that in order to obtain a 
satisfactory description of the data it is necessary to introduce a gluon coupling to gi. This 
is due to the fact that the SMC and E143 experiments measure gi at different values of 
for equal values of x. The two data sets can then be made consistent with each other only if 
there is a sufficient amount of perturbative evolution, and this can hardly be obtained in a 
LO computation^^) , or if schemes where the gluon decouples from gi are adopted. 

Inspection of the deuteron predictions (fig. 3b), however, clearly shows that more definite 
conclusions may be reached by combining the two sets of data: for instance the "minimal 
gluon" case seems to be excluded. This is a consequence of the fact that the deuteron is a 
more sensitive probe of the singlet and thus of the gluon, as explained in sect. 3.1. The results 
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Table 2: Next-toleading order determination^^ of polarized first moments. 



of a full NLO determination^) of gi (displayed in fig. 4) indeed confirm that all the parameters 
in eq. (4.8) may then be determined with reasonable accuracy. In particular, it is possible 
to disentangle the singlet and nonsinglet contributions, and thus determine independently 
the small x behavior of the various parton distributions. This turns out to be steep for the 
nonsinglet (Ag^g ~ x"^^ with q;ns = —0.7 ± 0.2), and valence-like for the singlet (with 

large error). The result is quite suggestive in view of the discussion in the previous section: 
it seems consistent with the expectation, based on the resummations of ref. ^'^\ that the 
polarized and unpolarized nonsinglets should behave in the same way, despite the prediction 
of Regge theory. However, the asymptotic small x behavior is only starting to set in at the 
smallest values of x covered by present-day data, which, in this respect, are to be compared 
with the NMC unpolarized data ^^^r a kinematic coverage comparable with that available 
at HERA for determining "^^^ the small x behavior of F2 would be required to determine 
precisely the small x behaviour of polarized parton distributions. 

Present-day data are however sufficient to determine the singlet first moments AE(1, Q^) 
and Agf(l,(5^) with reasonable accuracy (see table 2), assuming isospin [and SU(3)] to be 
exact. Even though the first moment of Fi at each determines only the linear combina- 
tion eq. (3.15) of AE(1,(5^) and A5f(l,Q^), their different scale dependence is sufficient to 
determine them independently: the large evolution effect already observed requires a large 
gluon component. Equation (3.15) then forces a large value of AE(1,(5^), consistent with 
the Zweig rule. Notice that the direct gluon contribution to gi is essential in order to get 
good agreement with the data. Thus current data support the proposal of ref. ^'^^ that the 
smallness of oq is due to a large gluon contribution to it. 

The first moment of the quark and gluon distribution then determine both the axial 
charge and the first moment of Fi. These determinations turn out to be significantly smaller, 
and with significantly larger error, than the corresponding determinations from the exper- 
imental collaborations (table 1): as discussed in sect. 3.1, a difference of the order of 10% 
in the value of F^ results in a variation of a factor of about 2 in the value of ao- We may 
actually understand in detail the origin of this discrepancy^°\ which is due to the approx- 
imate treatment of both scale dependence and extrapolation at small x used in refs.^~^). 
Firstly, the data points which provide the bulk of Fi are taken at a scale larger than the 
nominal average scale of either experiment: for example in refs.^'"^-' Qexp always larger 

than 20 GeV^ for x above 0.1, while (<5^) = 10 GeV^ for this experiment. Since in this region 
gi increases as the scale increases the substantial underestimate of evolution effects due to 
their approximate treatment in refs.^~^) leads to an overestimate of Fi. Furthermore, the flat 
small X extrapolation based on Regge theory overestimates the corresponding contribution, 
because it is obtained by extrapolating the smallest-x data points, which are taken at very 
low (~ 1 GeV^) but assumed to apply at larger values of where perturbative evolution, 
as explained in sect. 4.1, leads rapidly to negative values of gi and thus to a small (or even 
negative) contribution to Fi from the small x region. 
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The presence of large corrections due to perturbative evolution then inevitably implies 
sizable uncertainties related to insufficiently precise knowledge of the polarized gluon distribu- 
tion, which drives this evolution at least at the level of the first moment, and to the unknown 
higher order corrections (scheme ambiguities). These are indeed the dominant sources of the 
error on Fi and thus ao in table 2. It is interesting to notice that even though the values of 
table 2 are obtained by assuming isospin [i.e. the Bjorken sum rule eq. (2.8)] to be exact, the 
values of ao, AE(1) and Ag{l,Q'^) are essentially insensitive to SU(2) [or SU(3)] violations. 
This suggests that if the isotriplet first moment Fi were also determined by a global NLO 
analysis, rather than being fixed using the Bjorken sum rule, its value might not differ much 
from that found by combining values of Fi from table 1. Such an analysis would however be 
required in order to test the Bjorken sum rule in a reliable way. 



5. Outlook 



Polarized structure functions determined in deep-inelastic scattering are clearly the pre- 
ferred source of information on the polarized structure of the nucleon at high energy, since the 
availability of renormalization group methods and factorization theorems that characterizes 
fully inclusive processes allows a determination of physical observables which is free from the 
ambiguities due to our ignorance of the low-energy dynamics. In particular, due to the axial 
anomaly, polarized deep-inelastic scattering provides an ideal probe of the polarized gluon 
distribution and allows us to get a handle on quantities that have direct relevance for our 
understanding of the chiral structure of the QCD vacuum. 

The new generation of experiments^"^'' has allowed us to confirm these expectations, by 
providing us with surprisingly accurate information on polarized parton distributions. The 
analysis of these data, however, requires sophisticated theoretical and phenomenological tools. 
A future generation of experiments with a wider kinematic coverage in x and Q^, such as 
those that would be possible with a polarized proton beam at HERA^'^) , could lead to a full 
understanding of the small x behavior and scale dependence of polarized parton distributions, 
thereby testing our understanding of QCD in perturbation theory and beyond. 
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